Dynamics of defect formation 
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A dynamic symmetry-breaking transition with noise and inertia is analyzed. Exact solution of the 
linearized equation that describes the critical region allows precise calculation (exponent and prefac- 
tor) of the number of defects produced as a function of the rate of increase of the critical parameter. 
The procedure is valid in both the overdamped and underdamped limits. In one space dimension, 
we perform quantitative comparison with numerical simulations of the nonlinear nonautonomous 
stochastic partial differential equation and report on signatures of underdamped dynamics. 
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When a system that undergoes a symmetry-breaking 
transition is swept through its critical point the initial 
symmetry is broken and domains are formed. Because 
of critical slowing down it is not possible to sweep adi- 
abatically; the number of domains therefore depends on 
the rate of increase of the critical parameter. A new sce- 
nario for structure formation in the early universe and a 
proposal for its test in laboratory experiments resulted 
from the first understanding of the importance of this 
nonequilibrium effect jj| . Until recently, experimental j^[ 
results tended to support the proposed scenario, but pre- 
cise comparison was not possible because neither exper- 
iment nor theory was confident of more than exponents. 
The situation is now changing, with new experiments 
using quenches of liquid helium through the superfluid 
transition taking care to minimise vortex creation via 
flow processes ||. In this Letter we report new theoreti- 
cal results: precise expressions for the number of defects 
and quantitative agreement with numerical results. 

The phenomenon of a dynamic transition has been 
studied in the zero-dimensional case (pitchfork bifur- 
cation) in the context of lasers Q|- The time- 
dependence of the critical parameter produces a delay of 
the bifurcation given by •v/2/x] lne| where /i is the rate of 
increase of the parameter and e the magnitude of additive 
fluctuations. Theoretical studies on spatially-extended 
systems revealed a characteristic distance between kinks. 
The spatial structure formed during the sweep through 
critical point from the symmetric to broken-symmetry 
regime is frozen in by the nonlinearity when, sufficiently 
far into the symmetry-broken regime, the system at- 
tains a metastable state |H)- |l4| , |l6f . Analytical progress 
is possible because the critical region is well-described 
by an equation which, although stochastic and non- 
autonomous, is linear. Here we consider the influence 
of inertia; we derive the scalings and signatures of the 
overdamped and underdamped limits. 



cessive regimes in the evolution, as the critical parameter 
is increased. In the earliest regime, sufficiently far from 
the critical point, the evolution is quasi-adiabatic: the 
ensemble of field configurations is a small perturbation 
of that found for constant parameters [0,[l6) ■ In the sec- 
ond region, close to the critical point, the system can 
no longer react quickly enough to the time-dependence 
of the critical parameter Our treatment based on 
the equation of motion, however, passes seamlessly be- 
tween the first and second regions: in both, the field is 
everywhere small and precise calculation of the correla- 
tion function can be made from the linearized stochastic 
partial differential equation (SPDE) [0]l(|. We show 
that, for the purposes of calculating the number of kinks 
formed, the end of the second, nonequilibrium, region is 
the key. In the final region, the spatial structure con- 
sists of narrow kinks separating long regions where the 
field is close to one of the minima of the potential. The 
spatial structure is "frozen in" in the sense that the mo- 
tion, merging and occasional nucleation of kink-antikink 
pairs happens on a slower timescale than the process that 
formed them. The separation of timescales is especially 
marked at high damping and low temperature JL3|] . 

We shall consider the specific example of the stochas- 
tic process in one space dimension satisfying the following 
nonautonomous SPDE |T2|,f4|Jl5 1 



d?Y t {x) - Dd 2 x Y t {x) + jd t Y t (x) = 
= g(t)Y t (x) -Y t {xf + et]{x,t). 



(1) 



The theory of dynamic transitions identifies three sue- g 



The order parameter at time t and position x, denoted 
by Y t (x), is a real- valued random variable. The last term 
in ([!]) is space-time noise, delta-function correlated in 
space and time: (rj{x,i)rj(x' , i')) = S(x — x')S(t — t'). 
The fluctuation-dissipation relation is enforced by setting 
e 2 = 27fcs<9, where is the temperature. 

In our numerical simulations, the time-dependence 
of the critical parameter is: g(t) — [it, starting at 
-t < 0. The initial conditions are 



f 



Y to {x) = d t Y to (x)=0, t 



T 
(1 



(2) 



The simulations are performed on a domain [0, L] that 
contains many kinks, using periodic boundary conditions. 
Second order stochastic time-stepping [|17[ was used for 
the spatially discretized version of ([!]) |ll| . 

Typical time evolution is represented in Fig. [l], where 
each dot is the spacetime position of a zero crossing in 
one numerical realization. The system makes a transition 
from a regime with many zero crossings and typical val- 
ues of the field close to zero, to a regime with few zeros, 
corresponding to the positions of kinks, separating large 
regions where Y t {x) is close either to +y/g or to —y/g. 
The transition takes place at g > 0. For g > g, each zero 
of Yt (x) corresponds to a well-defined kink or antikink. 
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FIG. 1. Space-time evolution. The positions of crossings of 
zero are shown as a function of time for underdamped (a) and 
overdamped (b) dynamics. In (c) and (d) the corresponding 
numbers of crossings are shown. (Time increases upward.) 
The horizontal lines correspond to g — g. 

During the evolution preceding g = g, the cubic term 
in (Q) is small and the linearized equation is a good ap- 
proximation. It is illuminating to nondimensionalize the 
linearized equation: 



d 2 Y T (x 



where T 



vdlY T {x) 

= TYr(x) + efi 
fih, 



a = \-yp 



+ 2ad T Y T (x) 

-2T](x,t), 

3 , v = Dp~ 



(3) 



The dynamics can now be studied in terms of the char- 
acteristic time t = and nondimcnsional damping a. 
The field Yr(x) satisfying (^) with initial conditions (j^) 
is Gaussian with mean zero at all times. The correlation 
function, c(x) — (Yt(x)Yt(O)), changes its form and am- 
plitude with time. At any fixed time, there is the follow- 
ing relationship between c(x) and the number of zeros: 
if c'(0) = then the mean number of zero crossings is a 
finite number given by 



P{T) = n {- W 



(4) 



Analytical solution of (g) proceeds by separating into 
independent stochastic differential equations for each of 
the Fourier coefficients, t/y(fc) = J Q L Yr(x)e lk2 S' x dx, 
whose time evolution is given by the SDE 



for 



<9 T y T (fc) + 2adTU T (k) — 
= (T - K 2 )y T (k) + (2n)-hf,(T,k) 

k and where 



integer k and where k a = vk 2 {^) 2 
(fi(T,k)fi(T',k')) =6k 3 k>5(T-T'). Each 2 / T (fc) is Gaus- 
sian with mean zero ]19|. The variance grows exponen 



(5) 
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tially fast for T 

(y* T {k)y T (k)) - 
x exp (§(T- 
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ol z > 1: 
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$(T ,a,K) = e^ a3 / Ai 2 (S-K 

J To 



-K 2 + a 2 )e 2a »dS. (7) 



No approximations have been made thus far in the 
solution of (||). We now consider the implications of the 
physical picture presented above for the relative values of 
the parameters. Firstly, for there to be a quasi-adiabatic 
first regime in the evolution, we require a sufficiently slow 
sweep: /J, <C [ po[ . Secondly, we require a well-defined 
value of the order parameter, g = g, marking the end of 
the second part of the evolution, implying e 2 <C /i [flTl . 

We adopt the following definition: g = fi^T where T 
satisfies (Y?(x)) = 5g. (i.e., when 6 = 1, the first two 
terms on the RHS of (0) are equally important.) We thus 
evaluate g by solving (Y?(x)) — 6g where 



c(A(T))- 1 (27r)^(^(0)y T (0)) 



(8) 



and A(T) = -2v-§^ In ({y* T (k)y T (k))) . 

The correlation function is the Fourier transform of 
(y T (k)y T (k)} . It emerges from the sweep past g = 
with the form c(x) = c(0) exp(-a; 2 /2A 2 (f )) Q. The 
number of zeros present at g = g is thus 



I L 
p = — . 

7T A(T) 



(9) 



Our procedure is valid for arbitrary damping. We now 
examine the underdamped and the overdamped limits, 
defined by the parameter a = ^ 7/1^3 . The overdamped 
limit (studied in ]TT[| ) is recovered as a — > 00. 

The underdamped limit. When a — > 0, p is only loga- 
rithmically dependent on a. In this limit A(T) = 2v^ T', 
and the integral ([?]) has the asymptote <&(T , a, k 2 ) — > 
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and g satisfies 



? T-icxp(^ri), 
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Immediately before g — g, the number of zeros 
decreasing function of time, given by 



piT) 
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The number of zeros present at g for a — > is 
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is a 
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(13) 



In the overdamped limit the number of zeros is propor- 
tional to pi (l^Jl^]. We show that the latter scaling is 
obtained in the limit a — > oo. Here $(Tb, a, k) — » $2Q ;_ ^ 
H and A(T) 2 = 2i/T/a. Thus 



e 2 $ 2 i 1 T 2 
7T -U%vT)~? exp( ) (14) 



(y|(x)) 

and g satisfies 

g 2 = /i 7 ln(e-V(8^5 3 )^) 



(15) 



The number of zeros for overdamped slow passage is 



P(f) 



2tt 



1/4 



ln(e- 2 7(5(8Dg 3 )2) 



(16) 



Equations (|13|) and (|16|) are the main results of this 
Letter: the number of created defects scales with the 
sweep rate as p ~ p* for the underdamped case and 
p ~ /i3 for the overdamped regime. We have performed 
extensive quantitative comparison with numerical simu- 
lations of (^). Two examples are shown in Fig. ^. Our 
analytical predictions at instant g = g are very accurate. 
Although the exponents can be obtained from dimen- 
sional analysis pi[lq], logarithmic corrections produce 
small deviations in numerically-estimated exponents at 
finite damping. No evidence has been found for the re- 
gion of p ~ /i2 scaling, predicted in ]l5|j from an approx- 
imation that replaced (j^) by a first-order equation. 

Because the number of defects at g is typically much 
larger than the equilibrium density at temperature 
p2j , their number decreases after g as kink-antikink pairs 
annihilate (see Fig. [I]). The smaller the damping, the 
more rapidly this annihilation proceeds [ p"3"| . In Fig. ^ we 
have also plotted the number of zeros at g = r. While 
the number of zeros is reduced, the scaling with p seen 
a t 9 = 9 is preserved. 



The crossover between regimes, represented in 
Fig. pL takes place when the nondimensional damping 
a = 7py/i~3 ~ 1. At small damping the dependence of 
the number of defects on damping is only logarithmic, 
(e 2 oc 7). At large damping, p ~ 74. 
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FIG. 2. Density of zeros for the underdamped (a) and over- 
damped (b) cases. The lines are the analytical approxima- 
tions ( p^[ ) and (p^[). Averages from numerical simulations of 
([!]) are shown at g = g (squares) and g — r (circles), with 
errors of symbol size or smaller. {D = r = 1, = 5 X 10 -9 .) 
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FIG. 3. From the underdamped to the overdamped regime: 
density of zeros vs a for fixed p. Squares are numerical av- 
erages. The main plot shows the density of defects at the 
instant g and the analytical predictions for the underdamped 
( [l3|) and the overdamped ( llq) regimes. We measured g nu- 
merically as the instant when (Y?(x)^ = \g. In the inset we 
plot the value thus obtained vs ([H]) and (|l5|). (p — 10 -2 , 
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Apart from the scaling of the number of defects with 
/i, a different signature of underdamped dynamics can 
be seen in in Fig. |l| and Fig. [|: multiple "bounce back" 
of the number of zeros soon after g = g. The phe- 
nomenon has been reported in simulations of a sudden 
quench (p — 00) p3| ]. We propose the following interpre- 
tation. In a dynamic transition at low damping, domains 



3 



reach a minimum of a potential well with a finite veloc- 
ity and therefore oscillate about it for a time. Parts of 
some domains recross the crest of the instantaneous po- 
tential barrier during these oscillations. This yields an 
estimate of the frequency of the oscillations: \/2g, cor- 
responding to harmonic oscillations about the minimum. 
In our simulations, two well-defined bumps are typically 
seen in the number of defects vs time. From this we are 
able to measure the period, P of the oscillations in the 
number of zeros; despite the nonlinearity, it is very well 
approximated by P = 2n/^/2g (See Fig. 

The procedure carried out in this Letter for the real 
equation ([!]) can be applied to other equations exhibiting 
continuous transitions [fill and in more than one space 
dimension (l^jlj]. The scalings are not sensitive to the 
particular equations chosen, but they are sensitive to any 
breaking of the exact Y —> — Y symmetry in the equation 
of motion. 
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FIG. 4. "Bounce back": The number of zeros oscillates af- 
ter g = g for small damping. The graph of r vs time is for 
IX = 4 x 10~ 3 , 7 = 10~ 4 , e = 10" 6 and D = r = 1. The inset 
plots the period of the oscillations obtained at different values 
of jit, with 7 = 10 -4 . The straight line is P — 2n/^/2g with g 
given by (y_lj). 

In summary, we derive quantitative predictions for the 
number of defects formed in a symmetry-breaking tran- 
sition in one space dimension by analyzing the dynamics 
in the critical region, where the system is out of equilib- 
rium regardless of how slowly the critical parameter is 
changed. Prefactors are calculated, so no fitting neces- 
sary. Underdamped slow passage results in a defect den- 
sity proportional to /js and produces characteristic oscil- 
lations in the number of zeros. Experiments where liquid 
Helium is expanded through the Lambda Transition re 
now reaching the point where quantitative comparisons 
can be made. 

We are grateful for Angel Sanchez's comments on the 
manuscript. E. Moro thanks the CNLS for its hospitality. 
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